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Abstract 

The discrepancy between experimental and theoretical total electron-impact ionization cross- 
sections for a group of fluorides, CF X and NF X (x=l-3), is attributed to the inadequacies in 
previous theoretical models. Cross-sections calculated using a recently developed siBED model 
that takes into account the shielding of the long-range dipole potential between the scattering 
electron and target are in agreement with experiment. The present study also carefully re- 
analyzed the previously reported experimental data to account for the possibility of incomplete 
collection of fragment ions and the presence of ion-pair formation channels. For NF 3 , our 
experimental and theoretical cross-sections compare well with the total ionization cross-sections 
recently reported by Haaland et al. in the region below dication formation. 
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1. Introduction 

Electron-impact ionization of atoms and molecules is a fundamental process in any plasma. 
Large experimental databases for electron-impact ionization cross-sections for a number of 
molecules are now available [1]. However, measurements of reactive species such as radicals are 
known to be difficult. Essentially all available measured ionization cross sections for free radicals 
were obtained by the fast-neutral-beam technique which was pioneered by Freund and co- 
workers [2] and subsequently used extensively by Becker and co-workers [3], Theoretically, ab 
initio calculations present a challenging problem and currently are limited to atoms. For 
molecules, two physically based models, the binary-encounter Bethe (BEB) model of Kim and 
Rudd [4] and the semiclassical Deutsch-Mark model (DM model) [5] have provided efficient 
means for the determination of ionization cross-sections. The cross-sections determined using 
both models are in generally good agreement with experiment and with each other. However, for a 
group of fluorines, CF X and NF X (x— 1-3), discrepancies between 50% to 100% are found when 
compared with the experimental data of Tamovsky et al. [6-9], This is a surprising result in view 
of the good agreement obtained for many molecules. Also, the ionization cross-sections of NF 3 
from the recent measurement by Haaland et al. [10] using Fourier Transform Spectroscopy are in 
good agreement with the data of Tamovsky et al. [8] except in the region of dication formation. 
This result lessens the probability that experimental problems are responsible for the observed 
discrepancy. Up to the now, no adequate explanation has been proposed for the difference 
between theory and experiment. 

Both CF X and NF X are important gases in semiconductor processing. All three CF X are very 
reactive. They and their ions are responsible for most gas phase and surface reactions in 
fluorocarbon plasmas. Nitrogen trifluoride is used both in plasma etching and in thermal cleansing 
of semiconductors and liquid crystal display panels. In view of the importance of these fluorides 
in industrial applications, we have earned out a combined theoretical and experimental study of 
the electron-impact ionization of these molecules. This study should also improve our 
fundamental understanding of electron-molecule collisions. 

The group of molecules in this study shares two common features: they only have fluorine 
bonds and the non-fluorine atom has nonbonding valence electrons. While the charge distribution 
of a nonbonding electron usually centers near the atomic nucleus, in this case the strongly polar 
fluorine bond pulls the electron into the bonding region, resulting in a built-up of charge there. In 
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a theoretical model of electron collisions with this group of fluorides, the strong repulsive field in 
the bonding region must be taken into account in order to describe the collision process properlv . 

Recently one of us has developed an improved binary -encounter dipole (iBED) model [10] for 
electron-impact ionization of atoms and molecules. In this model, the long-range dipole 
interaction between the electron and target is represented by the dipole Bom cross-section, 
instead of the dipole Bethe cross-section as in the BED model. Furthermore, the dipole potential 
is properly shielded as the electron approaches the bonding region. When the optical oscillator 
strength (OOS) for the dipole transition is not known, a simplifed version of the iBED model 
(siBED) is introduced that uses an approximate OOS based on the f-sum rule. If the inadequate 
description of the dipole interaction potential is indeed the source for the discrepancies between 
theory and experiment, then the iBED/siBED model should be capable of reconciling the 
difference. 

A brief outline of the iBED/siBED model is presented in Sec. 2 together with the determination 
of the shielding and molecular parameters. The re-analysis of the experimental data is presented 
in Sec. 3. The experimental and theoretical cross-sections are compared in Sec. 4 and Sec. 5 
concludes our results. 


2. Theoretical treatment 

Consider the spherical harmonics expansion of the Coulomb potential 
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The dipole term is given by 
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When the scattering electron is outside the molecular charge cloud, the dipole potential is 
described by 
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with D u one component of the transition dipole moment. When the scattering electron is inside 
the molecular charge cloud, the interaction potential changes to 

IVO'XW- <4> 


where B^frJ is given by 
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Here p po (r,Qj is the transition charge density of the target. Eqs. (3) and (4) show that the dipole 
potential changes character when the scattering electron penetrates inside the molecular charge 
cloud. The potential in Eq. (3) is the long-range dipole potential, and the short-range potential m 
Eq. (4) can be considered a shielding potential because it prevents the long-range potential in Eq. 
(3) from diverging at r 0 = 0. Note that the dipole Bethe cross-section used in the BED model [3] 
to describe dipole interaction is a high-energy approximation. At the limit of very high energies, 
only glancing collisions are important. Thus the dipole Bethe cross-section only accounts for the 
long-range dipole interaction in Eq. (3), but not the short-range shielding of the dipole potential in 
Eq. (4). 

For the ionization of an electron out of the orbital o in the target to the ion state p, the iBED 
cross-section is given by 
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Here O p( f B is the dipole Bom cross-section and <y po BE is a modified Mott cross-section with the 

incident electron energy replaced by the average energy from the binary-encounter model. The 
Mott cross-section is a generalization of the Rutherford cross-section for Coulomb scattering by 
taking into account exchange. The direct term in the Mott cross-section is associated with the 
1=0 term in Eq. (1). a po BE is referred to as the binary-encounter cross-section by Kim and Rudd 
[4] but is different from the binary-encounter cross-section of Vriens [1 1], It is given by 
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Here N a is the occupation number of the orbital, cc 0 is twice the magnitude of its binding energy, 
twice its kinetic energy, and k 0 the momentum of the incoming electron. 

The dipole Bom cross-section is represented by a three-term polynomial derived from a series 
representation for the generalized oscillator strength (GOS) of electron impact ionization. The 
series representation for the GOS is the ionization analog of the Lassettre series for bound-bound 
transitions [12]. This representation not only describes the long-range dipole interaction of the 
electron and target, but also the shielding of the dipole field as the scattering electron comes inside 
the bonding region. The use of the dipole Bom cross-section also remo'ves the empirical energy- 
scaling factor in the BED model. The following expression for the dipole Bom cross-section is 
used in the iBED model, 
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Here E p and k p are the energy and momentum of the ejected electron, W po the excitation energy, K 
the momentum transfer of the scattering electron, and df po (o> (Ep)/dE p the optical oscillator strength 
(OOS) for the corresponding photoionization process. The parameters dj and <7; are related to the 
shielding of the dipole potential, and the variable t is given by 
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The total single ionization cross-section of a molecule is given by 


a =YYe' BED . 0 °) 
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The calculation of iBED cross-sections requires the OOS for the corresponding photoionization 
process, the binding energy and the kinetic energy of the electron being ejected, and the shielding 
parameters dj and d 2 . In the case when the OOS is not available, an approximation expression has 
been proposed based on the f - sum rule [10]. 
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In the simplified version of the iBED model (siBED), Eq. (1 1) is used to represent the OOS. The 
siBED model is used for all CF X and NF X calculations reported in this paper. 

2. 1 Choice of d ) and d 2 

In principle, d, and d 2 can be calculated using the interaction potential in Eq. (4), with BfrJ 
determined using Eq. (5) and the transition density p po deduced from target wave functions. 
However, for molecular wave functions represented by a Gaussian basis set, the partial integral in 
Eq. (5) is not in closed form. In Ref. 9, a different approach is used to determine d) and d 2 . For 
N 2 and CH 4 , the two parameters were determined empirically based on the optimal 
representation of the experimental data of Straub et al. [13, 14] and Rapp and Englander-Golden 
[15] using the iBED model [10]. (Note that a recent recalibration of the apparatus by the Rice 
group results in data that differ slightly from those in the original publication [16]). For both N 2 
[10] and CO [17], the optimal choice of the two parameters are d, = -2.0, d 2 = 0.4. The 
parameter d, is found to be important to in determining the shape and magmtude of the cross- 



sections near the peak of the cross-section versus energy curve. The negative value denotes a 
repulsive correction to the purely attractive long-range dipole potential as the scattering electron 
enters the molecular charge cloud. The parameter d 2 does not play any significant role until the 
incident electron energy is above 100 eV. It provides a small attractive correction to d, as the 

scattering electron moves close to the nuclear centers. 

In the present study, we determine dj and d 2 by comparing the molecular property f e for CF 

and NF with CO. The property f e is defined by 
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Note that the/ e closely resembles the expression BJrJ in Eq. (5) that determines the short-range 


shielding potential. While B/rJ uses the transition density p po and the integration over r, is over 
a limited range, f e is defined by the density p DO of the initial target state and the r, integration 
covers the full space. How f e changes from CO to the fluorides may serve as an indicator how 
B^rJ, and consequently d h should be modified. Table I compares the values of f e for the valence 
orbitals of CF and CO, and for NF and CO. The orbitals are determined using complete-active- 
space self-consistent-field (CASSCF) calculations. These calculations used the aug-cc-pVQZ 
basis of Gaussian functions [18]. The active space used in the CASSCF calculations is (6, 3, 3,0) 
for CF and NF and (6, 2, 2,0) for CO. To eliminate variations in f e due to the difference in 
geometry, the experimental geometry of CF is used in the calculation for the CF/CO comparison 
and the NF/CO comparison uses the experimental geometry of NF. Table 1 shows that the 
orbital f e 's for CF and NF, with one exception, are close to a factor of 2 larger than the 
corresponding f e ’ s for CO. Based on this result, and similar results from the comparison of the 
orbital s of CF 2 and NF 2 with C0 2 , and CF 3 and NF 3 with BF 3 , we determine that the increase 
shielding in this group of fluorides can be represented by scaling the d, parameter by a factor of 
2.4. Thus d, = -4.8 is used in our fluoride calculations. 

There is little guidance available in choosing the parameter d 2 . However, the role of d 2 is 
insignificant until the incident electron energy is above 100 eV. We slightly modify the d 2 of CO 
to reflect the diminished attractive potential in the inner region of the fluorides and use d 2 = 0.3. 
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2.2 Calculation of molecular parameters 

The binding energies and kinetic energies of the target electrons have been determined usmg 
quantum chemistry calculations. All calculations use the experimental geometry of the neutral 
molecules. The calculations for CF, CF 2 , NF, and NF 2 , employ the aug-cc-pVQZ basis of 
Gaussian functions and the aug-cc-pVTZ basis for CF 3 and NF 3 . The kinetic energies of the 
target electrons, k 0 2 /2, are directly calculated using Hartree-Fock functions. On the other hand, 
the binding energies -a Q 2 /2 are sensitive to correlation effects. For valence electrons, they are 

determined by taking the energy difference between the neutral target state and the ion. Smce the 
calculations are done at the equilibrium geometry of the target, they are just the vertical ionization 
potential (VIP) associated with the ion states. The RCCSD(T) method (symmetry restricted 
couple-cluster singles and doubles with perturbation treatment of triples) [19,20] is used m the 
calculation of the lowest ion state of each symmetry. The RCCSD(T) method is size consistent 
and thus suited for the calculation of energy difference between systems with different number ot 
electrons. However, RCCSD(T) can only be used to determine the lowest ion state of each 
symmetry. The energies of the second and higher ion states of a symmetry are determined using 
CASSCF calculations and the result is scaled by the difference between the CASSCF and 
RCCSD(T) result for the lowest state. The binding energies of the inner orbitals are determined 
from Hartree-Fock functions using Koopmans theorem. All quantum chemistry calculations have 
carried out using MOLPRO [21]. Table 2 presents the binding energies and kinetic energies for 
the fluorides. The data for the core orbitals are not listed because their contributions to the 
ionization cross-section are negligible in the energy range considered in this study. For CF 2 , the 
binding energies of the 4a,, 5a,, and 3b 2 electrons are taken from the compiled ionization potential 

data of Levin and Lias [22]. 


3. Analysis of experimental data 

We carefully re-analyzed the originally reported experimental cross section data for all six 
free radicals in an effort to ensure that systematic effects that may influence the proper 
determination of the total single ionization cross sections are fully accounted for to the maximum 
extent possible. The total single ionization cross section for each free radical is obtained as the 
sum of all measured partial ionization cross sections. The sum of all partial ionization cross 
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sections for a given target will only yield the “correct” total single ionization cross section, if all 
channels leading to the formation of one singly charged positive ion are properly taken into 
account. This requires estimates for contributions to the total single ionization cross section 
from partial ionization channels whose cross sections were not or not accurately measured 
because of weak signals and/or significant excess kinetic energy that precludes a 100% collection 
of a specific fragment ion. Another issue is the presence of (positive) ion - (positive) ion pair 
formation processes whose contributions to the total single ionization cross section have to be 
discounted which requires that their respective absolute cross sections and energy dependences 
have to be determined or estimated. In the course of the present work we re-analyzed all original 
data files relating to the measurements of the CF X and NF X partial ionization cross sections in an 
effort to obtain the most reliable experimental total single electron impact ionization cross 
sections for these six free radicals. 

4. Results 

Figure 1 presents the total single ionization cross-sections calculated using the si BED model 
and the revised experimental values for CF, CF 2 , and CF 3 . Also presented are the theoretical 
cross-sections using the BEB model and the DM model. The corresponding data for NF, NF 2 , 
and NF 3 are presented in Fig. 2. For NF 3 the data from the recent experiment of Haaland et al. 
[23] are also presented. Note that all BEB cross-sections have been calculated using the molecular 
parameters determined in this study. The cross sections therefore differ slightly from the BEB 
calculations reported by Kim and Irikura [24]. For all six molecules, the siBED cross-sections and 
the revised experimental cross-sections agree to within experimental error. For CF 2 , CF 3 , NF 2 and 
NF 3 , the agreement is very good and the two sets of cross-section curves almost lie on top of 
each other. The agreement is less satisfactory for CF and NF where the siBED cross-sections are 
consistently larger than the experimental values. Flowever, as pointed out in the original papers 
by Tarnovsky et al. [7,9], the atomic fragment ions resulting from the dissociative ionization of 
these two diatomic radicals are formed with significant excess kinetic energy. It is, therefore, 
difficult to quantify the collection efficiency of the fragment ions in these cases and the reported 
experimental ionization cross sections for NF and CF should be considered a lower limit [7,9]. 
Our cross-sections for NF 3 also agree well with the data of Haaland et al. [23] at 100 eV and 
below. It should be noted that the measurements of Haaland et al. includes the contributions from 
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all ionization processes, among them dication formation. For electron energies above the 
threshold of dication formation, the cross-section data of Haaland et al. should be larger than the 
corresponding total single ionization cross-sections. Thus their cross-sections and the cross- 
sections presented in this paper are no longer directly comparable at those energies. Also, as seen 
in Figs. 1 and 2, calculations using the BEB and DM models consistently overestimate the cross- 
sections for this group of fluorides. 


5. Conclusions 

The agreement between experiment and the siBED model suggests that the long-standing 
discrepancies between theory and experiment in the total single ionization cross-sections for CF X 
and NF X (x=l-3) are due to an inadequate description of the electron interaction potential in 
previous theoretical models. 
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Table 1. Comparison of the property f (au) for the valence orbitals of the CF/CO and for NF/CO 
pairs. Calculations are carried out at the experimental geometry of CF and NF, respectively. 


Orbital 

CF 

CO 

fe 

NF 

CO 

3a 

0.9665 

0.4738 

0.7230 

0.4947 

4a 

0.5595 

-0.0974 

-0.3091 

-0.1285 

5a 

-0.0577 

0.4364 

0.6524 

0.4383 

\k 

0.8296 

0.4152 

0.7372 

0.4070 
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Table 2. Binding energies (- a 0 2 /2 ) and kinetic energies (k 0 2 /2) used in the siBED calculations. 
Both quantities are in eV. 


Orbital 

CF (X 2 n) 

NF (X 3 A) 


-cto/2 

k 0 2 /2 

-a 2 /2 

k 0 2 /2 

3 a 

46.804 3 

101.378 

46.1 12 a 

98.804 

4a 

23.304 b 

80.062 

24.41 8 b 

78.282 

5a 

13,985 c 

48.144 

16.624° 

70.199 

1 K 

18.41 l c 

81.922 

17.913° 

84.018 

2n 

9.475° 

45.231 

12.651° 

60.703 

Orbital 

cf 2 

(X'A,) 

nf 2 

(X 2 B,) 


-cc 0 2 /2 

k 0 2 /2 

-a 2 /2 

k 0 2 /2 

3ai 

47.674 a 

95.446 

47.073 3 

92.903 

4ai 

24.000 d 

82.088 

29.489 a 

83.938 

5aj 

19.200 d 

75.808 

18.749 b 

74.653 

6a] 

12.244° 

55.719 

14.716° 

77.847 

lb] 

19.241° 

74.534 

18.380° 

79.830 

2b, 



12.195° 

68.900 

2b 2 

45.446 a 

106.832 

44.668 a 

106.706 

3b 2 

22.200 d 

81.676 

20.906 b 

78.274 

4b 2 

16.665° 

95.114 

16.273° 

93.850 

la 2 

17.298° 

88.008 

16.508° 

89.544 


Orbital 


CF 3 (X 2 A,) 


NF 3 (X'A,) 



- a ,;/2 

k 0 2 /2 

- cl , 1 ; /2 

k 0 2 /2 

3a, 

47.552 a 

93.051 

48.828 a 

87.299 

4a, 

26.400 a 

82.590 

30.523 a 

89.586 

5a, 

21.723 a 

75.316 

23.01 8 a 

69.615 

6a, 

1 1.056 c 

59.408 

1 3.73 l c 

81.082 

la. 

16.363 c 

94.458 

16.1 80 c 

96.263 

2e 

45.250 a 

105.136 

45. 1 27 a 

106.502 

3e 

22.983 a 

76.776 

23. 190 a 

75.853 

4e 

18.224 b 

85.933 

18375 b 

86.623 

5e 

17.590 c 

90.592 

16.443 b 

93.111 


a Hartree-Fock calculation based on Koopman’s theorem. 
b CASSCF calculation corrected by RCCSD(T) results. 
C RCCSD(T) calculation. 

Compilation of Levin and Lias [22]. 


Figure captions 

Figure 1. (a) Total single ionization cross-section of (a) CF, (b) CF 2 and (c) CF 3 calculated usmg 
the siBED model and the revised experimental data. Also presented are theoretical cross sections 
calculated using the BEB model of Kim and Rudd [4] and the DM model [5], 

Figure 2. (a) Total single ionization cross-section of (a) NF, (b) NF 2 and (c) NF 3 calculated using 
the siBED model and the revised experimental data. Also presented are theoretical cross sections 
calculated using the BEB model of Kim and Rudd [4] and the DM model [5], For NF 3 the cross- 
section data form the measurements of Haaland et al. [23] also are presented 















